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Abstract 

We show how to couple two free-fermion chains so that the excitations consist of Cooper 
pairs with zero energy, and free particles obeying (mutual) exclusion statistics. This behavior 
is reminiscent of anyonic superconductivity, and of a ferromagnetic version of the Haldane- 
Shastry spin chain, although here the interactions are local. We solve this model using the 
nested Bethe ansatz, and find all the eigenstates; the Cooper pairs correspond to exact-string 
or "0/0" solutions of the Bethe equations. We show how the model possesses an infinite- 
dimensional symmetry algebra, which is a supersymmetric version of the Yangian symmetry 
algebra for the Haldane-Shastry model. 

1 Introduction 

The Bethe ansatz is one of the few reliable methods physicists have for probing properties of 
systems with strong correlations. For the systems to which it is applicable, it gives a way of 
doing exact computations. 

One of the most valuable insights it gives is in the quasiparticle spectrum of one-dimensional 
quantum systems. A number of profound but highly non-obvious properties have been discovered 
using the Bethe ansatz. For example, Bethe's original paper gives the tools necessary to show 
that the quasiparticles in the one-dimensional antiferromagnetic Heisenberg model are gapless 
and have a linear dispersion relation [TJ. It took more than half a century and some novel non- 
perturbative physics to understand that these excitations have spin 1/2, and to explain this in 
terms of field theory [21 [3]. After the Bethe ansatz was generalized to systems with multiple 
types of particles via the invention of the nested Bethe ansatz [3], it was subsequently shown how 
the spectrum of the one-dimensional Hubbard model is spin-charge separated. Even though the 
Hubbard model is comprised of electrons with both charge and spin, the quasiparticles in one 
dimension have either spin or charge, but not both |5]. It was again quite some time before this 
result was understood in depth from the field-theory point-of-view. 

In this paper we use the Bethe ansatz to find and describe a strongly-correlated one-dimensional 
system with two kinds of interesting excitations. The system consists of two coupled free-fermion 
chains; the excitations are Cooper pairs and exclusons. 



Cooper pairs of course are familiar from studies of superconductivity and fermionic superflu- 
idity. They are comprised of two fermions bound (in momentum space) into a state which has 
zero energy. The Cooper pairs here arise in a fashion very reminiscent of anyonic superconduc- 
tivity [6]. Anyons are quasiparticles in two dimensions which have statistics generalizing those 
of fermions and bosons: they pick up a phase under exchange, instead of a sign. One kind of 
interaction between our fermions is a one-dimensional analog: a fermion on one chain picks up a 
phase ±i when it moves past a fermion on the other chain. 

Exclusons are also a generalization of the idea of fermions and bosons [7]. However, as 
opposed to anyons, they can occur in one (or any) dimension. The usual notion of statistics is not 
applicable in one dimension because particles can not be exchanged without coming close to each 
other. Nevertheless, in any dimension fermions can be defined as particles which are allowed at 
most one to a level, and bosons those where any number is allowed. Exclusion statistics describe 
more general rules for how particles fill energy levels. In our model, we will show that there 
are two types of exclusons, associated with each of the two chains. Each of the two types has 
their own (identical) set of energy levels, which with appropriate boundary conditions are those 
of free particles. Amongst themselves, each type behaves like a fermion, allowing at most one 
quasiparticle per level. However, when say an excluson of type 1 is in a given energy level, not 
only does it forbid another of type 1 from having that energy, but it also forbids a particle of 
type 2 from having that energy as well. 

Exclusons in one dimension are well known to occur as quasiparticle excitations of the 
Haldane-Shastry chain, which is a variant of the Heisenberg model with long-range interactions 
[8]. In both our model and the Haldane-Shastry chain [9], a nice property is that the interaction 
between exclusons is purely statistical. In other words, the energy levels are those of free par- 
ticles: the only interactions are in the rules for filling these levels. It is important to note that 
it is the quasiparticles which have only exclusion interactions: the original degrees of freedom in 
both cases are strongly interacting. One very interesting distinction between our model and the 
Haldane-Shastry chain is that the interactions here are nearest-neighbor. 

Even after 75 years, computations using the Bethe ansatz are usually quite technical, and this 
one is no exception. However, we find it remarkable that here the technical complications result 
in elegant physics. For example, to obtain our results we impose open boundary conditions, which 
require much more effort to deal with in the Bethe ansatz than periodic ones do. Nevertheless, 
not only does the Bethe ansatz still work here, but the fact eigenstates are standing waves ends 
up making it possible to find them all in closed form. This is how we manage to show that this 
strongly-interacting system has quasiparticle excitations which are free exclusons and Cooper 
pairs. 

The Cooper pairs themselves arise as a consequence of resolving another technical complica- 
tion: they turn out to be "exact-string" solutions of the Bethe equations. The Bethe equations are 
coupled polynomial equations whose solutions determine the allowed momenta in the wavefunc- 
tion. String solutions are those for which some of these polynomials vanish in the thermodynamic 
limit; they occur even in the Heisenberg model discussed in Bethe's original paper. In an exact- 
string solution, they vanish even for a finite number of sites. In the way the Bethe equations are 
usually written, this is a singular "0/0" limit: both the numerator and the denominator of a ratio 
are exactly zero. However, it is explained in detail in |10| how to define the Bethe equations so 
that exact-string solutions are not singular. Including such solutions is necessary to ensure that 
all eigenstates of the Hamiltonian follow from the Bethe ansatz. 

Moreover, the presence of exact-string solutions in the XXZ spin chain described in [10] is 
closely related to the appearance of a very interesting extended-symmetry algebra there |llj . 
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We find an analogous algebra for our model as well. Because our model turns out to have a 
supersymmetry, this extended symmetry algebra is an extension of the supersymmetry algebra. 
A similar symmetry algebra called the Yangian appears in the Haldane-Shastry chain |12j . 

In section [2j we present our model and describe in detail the results of this paper. In section 
El we derive these results using the Bethe ansatz. In section [J we discuss the extended symmetry 
algebra, and its connection to supersymmetry. We present our conclusions in section [5j 



2 The model and the results 



We consider two free-fermion chains, each with L sites. It is convenient to label the sites on 
one chain by odd integers 1, 3, . . . , 2L — 1 and on the other by even integers 2,4,..., 2L. The 
spinless fermions on the two chains are created by the operators cj which anticommute except 

for {cfc,cj} = 5jk- The fermions on a given chain do not interact among themselves, except of 
course for the requirement that two fermions cannot occupy the same site. The number operator 
n,j = CjCj is 1 if the site j is occupied and if it is unoccupied. We allow hopping along each 
chain, but not from one chain to the other, so the fermion- number operators 

L L 

k=l k=l 

are each conserved. A state with fermion numbers /i and f 2 on the chains is of the form 

v (ii, • • • , Jh b'/i+i , • • ■ i/i+/ 2 ) 4i • • • 4/1+/2 1°) ' W 

where |0) is the empty state. Here and for the remainder of the paper, j a is odd if a < /i and 
even if a > fx. 

Let us first consider two decoupled free-fermion chains with open boundary conditions. We 
include hopping, a specific chemical potential, and a boundary magnetic field, so that the Hamil- 
tonian is 

2L-1 

H f = 2F X + 2F 2 +Y, [ c ]+i c J-i + c ]-iw] + H b ■ ( 2 ) 
i=2 

The specific boundary magnetic field we choose gives 

H b = -ni - n 2 L ■ (3) 

so that it lowers the energy when a particle is on one end of each chain. 

An elementary computation gives the eigenstates of Hf to be of the form ([1]) with 

fl+f'2 

^free = II sin 0^) ( 4 ) 

a=l 

where 

Pa = TTla (5) 

with m a an integer from 1 . . . L. Even though the open boundary conditions mean that the 
system is not translation-invariant and there is no overall momentum conservation, we still refer 



3 



to the p a as the momenta. The shift in the denominator of © is due to the boundary magnetic 
field. The energy of this eigenstate is 

h+h 

E= (2 + 2cos(2p a )). (6) 

0=1 

To impose Fermi statistics and forbid double occupancy, one sums over these eigenstates as in 
(|12p below. One can have m a = only if a = b or if a and b label particles on different chains, 
i.e. a < fi and b > f±. For a given f\ and fi there are therefore (^)(^) different choices of 
momentum. This is the number of states, so (j3J) is a complete set of states. 

We now couple the two chains in a zig-zag fashion, so that each site j is adjacent to the 
sites j ± 1 on the other chain. There are two kinds of interaction between the chains. The first 
is simply an attractive interaction between nearest-neighbor fermions on different chains. The 
second is a phase factor dbi picked up when a fermion on one chain hops past one on the other, 
which can be thought of as resulting from flux tubes attached to the particles. The Hamiltonian 
is then Hf + H c , where 

22,-1 22,-1 

H c = - y~] [(1 - i)cj +1 rijCj-i + (1 + i)ct_ 1 n J Cj+i - 2 ^ . (7) 

i=2 j=l 

The specific values of the couplings in ([2|) and ([7]) are chosen because they result in a model with 
a great deal of symmetry, which we will discuss in depth in this paper. Some of the symmetries 
(e.g. integrability) persist under certain deformations of the couplings, but we will confine our 
studies to this specific Hamiltonian. 

The model with Hamiltonian Hf + H c is solvable using the nested Bethe ansatz, and in the 
next section we present the solution in detail. The Bethe equations can be solved explicitly, and 
the resulting physics is quite elegant. We therefore discuss the key physical results here. 

We find that in this strongly-coupled theory, the energy levels are still those of free fermions, 
i.e. they are given by ([6]) with momenta quantized in the free- fermion form ([5]). However, the 
theory is not the same as that of the decoupled chains: the strong interactions have several 
interesting effects. The wavefunctions can still be written in a form quite similar to (Jl|), so that 
each state can be characterized by a set of momenta p a with a = 1, . . . , fx + f-i- The energy is 
still given by ([6]). However, not all momenta obey ([5]) or are even real. A pair of particles, one 
from each chain, can form a Cooper pair. Namely, two of the momenta (say p\ and P2) obey 

cos 2 (pi) = -cos 2 (p 2 )- (8) 

This existence of Cooper pairs does not contradict the first sentence of this paragraph, because 
each pair obeying (|8|) has no net contribution to the energy. If f\ = fy, then one can pair all 
the particles and obtain a state with energy zero. We show on general grounds in section 0] that 
E > 0, so a state comprised entirely of Cooper pairs is a ground state. 

Any momenta not part of a Cooper pair must obey ©, so the energy levels are indeed those 
of free fermions. The degeneracy of each level, however, is not the same as that for the decoupled 
chains. One can add Cooper pairs without changing the energy, which obviously will change the 
degeneracies. Even in a state with no Cooper pairs, there is a crucial difference. In both coupled 
and decoupled chains, the numbers f± and / 2 of fermions on each chain are conserved. When 
the chains are decoupled, two momenta can be the same, if they are associated with particles on 
different chains. When the chains are coupled via H c , we find that no two momenta can be equal. 
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In fact, the individual momenta cannot be associated with one chain or the other like they can 
in the decoupled case: for C Cooper pairs, the wavefunction is labeled by /i + fa — 2C distinct 
momenta obeying ([5]) , and C pairs of momenta obeying (JBJ) . 

Thus we can interpret the states with non-zero energy in terms of quasiparticles with momen- 
tum @. We call these quasiparticles exclusons for the following reason. Consider fixed numbers 
fi and /2 of fermions on each chain, and also a fixed number of Cooper pairs C. There are two 
types of exclusons: Aq = f\ — C of type 1, and N 2 = /a — C of type 2. The energy of this state 
is given by 

N!+N 2 

E= 4cos 2 (m a vr/(2L + 1)) (9) 

a=l 

with the m a G {1, . . . L} distinct integers. We show in section [3] that the number of states with 
this energy is 

*=(^)( £ -r*) ™ 

We have checked this formula for the degeneracy directly by numerically diagonalizing the Hamil- 
tonian (which of course is how we originally discovered it). 

In the quasiparticle interpretation of (I9)10|) . the first binomial arises from the exclusons, the 
second from the Cooper pairs. The iVj exclusons of each type are indistinguishable from each 
other, so the number of ways of assigning Aq + N 2 momenta to them which result in distinct 
wavefunctions is ( Nl ^ 2 ) ■ The interpretation of the second piece is that the Cooper pairs can be 
thought of indistinguishable quasiparticles like the exclusons, and there is one possible Cooper 
pair for each allowed momentum. Moreover, the presence of an excluson with a given momentum 
prohibits a Cooper pair from occupying the corresponding level. 

The exclusons, the quasiparticles of non-zero energy, obey mutual exclusion statistics [7J. 
Even though in our counting we have treated exclusons on chain 1 and on chain 2 as distinct 
species, only one can fill a given momentum level. This is thus a generalized exclusion principle: 
the presence of one kind of quasiparticle affects how others fill levels. A mutual exclusion principle 
is typical for systems solvable by the nested Bethe ansatz; what is novel here is that except for 
this mutual exclusion, the quasiparticles are otherwise blind to each other's presence. Namely, 
the energy levels are unaffected by how many quasiparticles are present: the only effects are on 
how the levels fill. 

This mutual exclusion principle involves the Cooper pairs as well. Fermions of the same type 
must still of course obey Pauli exclusion. Thus the presence of exclusons restricts the allowed 
momenta of the Cooper pairs, and reduces the size of the space of states for the Cooper pairs as 
well. 

One important issue to note is that the momentum of one of the particles in a Cooper pair 
is arbitrary, as long as it is not the same as any of the exclusons. The arbitrariness is due 
to the degeneracy: any linear combination of the oIe states with fixed excluson momenta is an 
eigenstate of the Hamiltonian. This arbitrariness in choosing certain momenta in the Bethe 
ansatz is discussed in depth in [10] , and we will return to this issue in section [3) The counting 
rules become clearer if we choose any arbitrary momenta to be one of the values obeying ([5]) 
which are not already present. The momentum of the second particle in each Cooper pair is then 
fixed by ([8]). The second factor in the degeneracy (I10p then simply corresponds to the number 
of different such choices. 

The factor oIe looks somewhat complicated. However, it ends up yielding an extremely simple 
and elegant formula for the grand canonical partition function. The allowed momentum values 
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© are independent of how many quasiparticles are present, so we can treat the grand canonical 
partition function as a product over momentum levels in the usual free-particle fashion. The 
reason this is possible here is that the way to fix the arbitrariness in the Cooper-pair momenta 
uniquely is to make it obey the same rules as the exclusons do. In other words, there are four 
possibilities for each of these momentum levels: it is empty, occupied by an excluson on chain 1 or 
by an excluson on chain 2, or occupied by a Cooper pair. The first and last of these possibilities 
have energy 0, whereas an excluson in the ath momentum level has energy 4cos 2 (p a ). If Ai and 
A2 are the fugacities for each fermion on chains 1 and 2 (not including the chemical potentials of 
—2 already in ([2j)), the grand canonical partition function at temperature T is then 



Y[ [l + A c + (Ax + A 2 ) e - 4cos2 (W(2L+i))/T (n) 

771=1 



where Ac* = A1A2. In the limit T — > 00, Z simply counts the number of allowed states. This limit 
indeed yields the correct Z — > A L = 2 2L for two chains of length L (each site is either occupied 
or unoccupied). From (jllH one can compute any thermodynamic quantity, such as the specific 
heat, or the density of Cooper pairs. 

One can easily recover the degeneracies in (I10p by expanding out the product in Z, i.e. 



^=0^1=0 C=0 {m a } 

where M = Ni + i\T 2 is the total number of exclusons, E is given in Q and the sum ^{m a } ^ s 
over sets of integers obeying 1 < mi < m-2 < • • • < mj^ < L. Note that the maximum number of 
exclusons is L, so if f\ + / 2 > L, then C > f\ + f2 — L. 

One can find a Hamiltonian with the same grand canonical partition function by considering 
spin- 1/2 fermions with long-range interactions. Having two fermions of different spin and the 
same momentum corresponds to a Cooper pair. To make the Cooper pair have energy zero, 
the Hamiltonian must include a binding energy for fermions of the same momentum. Such 
a Hamiltonian is simple to write in momentum space, but when Fourier-transforming back to 
position space, one obtains a complicated long-range interaction. 

Before proceeding to the solution by Bethe ansatz we comment once again on the fractional 
statistics aspects of the model. We already observed that particles on a given chain act as 
fermions, while moving a fermion on one chain past one on the other involves phase factors of ±i. 
This then suggests a statistics matrix (in the sense of Haldane [7]) G = {gij} with gyy = 522 = 1, 
512 = 521 = 1/2 for these 'bare' particles. For the exclusons, which fully exclude one another 
from a given level, the natural assignment would be g\\ = 522 = 1> 512 = 521 = 1- We observe 
though that the 1-level partition sum that is associated with this choice of G (using results of 
|13l I14j) is different from the simple result that we find here. 



3 The solution using the nested Bethe ansatz 

In this section we derive all the claims made in the last section by using the Bethe ansatz. 

When all fermions are on the same chain, there are no interactions save forbidding double- 
occupancy. Exact eigenstates are given in iQ) above. To forbid double occupancy and anti- 
symmetrize them appropriately, we introduce a permutation P = (PI, P2, . . . , Pf) of the integers 
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(1,2,..., /), where / is the total number of fermions. Then we have 



/ 

= Et- 1 )^' II S ' m (PPaja) ■ (12) 
P a=l 

where |P| is the order of the permutation, and p a satisfies ([5]). 

The Bethe ansatz is a very natural generalization of (|12p . In situations like ours where there 
is more than one species of particle, imposing Fermi statistics does not automatically fix the 
relative coefficients between different orderings. One must therefore use the nested Bethe ansatz 
[4J. It is thus necessary to label the orderings by another permutation Q of (1, 2, . . . /), so that 
<pQ is the part of the wavefunction where jq\ < jQ2 < ■ ■ ■ < jqf- It is simplest to first work with 
periodic instead of open boundary conditions. This amounts to changing the sums in Hf in ^ 
and H c in ([7|) to run from 1 to L, interpreting positions mod 2L, and setting = 0. The Bethe 
ansatz for the eigenstate with periodic boundary conditions is then 




<pQ = E A Q P exp [i E PpjQa • (13) 



For open boundary conditions, the eigenstates consist of sums over these cp® with ±p a , i.e. are 
standing waves. Note that the momenta are permuted over all / = f% + fi fermions. 

The miracle of the Bethe ansatz in general is that in some situations, the p a and the A% can 
be found which result in an eigenstate of the Hamiltonian. The corresponding energy is exactly 
that given in ([6]). The p a are given as solutions of what are called the Bethe equations. The 
additional miracle here is that for open boundary conditions with ([3]), the Bethe equations can 
be solved explicitly, yielding the results described in section [2J 

Let us first discuss the case of two fermions. If both are on the same chain, the only constraint 
is no double occupancy, i.e. <p®(j,j) = for all Q. This means that 

A% = -Ag for h = 2, h = 0, (14) 

for both Q = (12), (21). The solution is much more interesting when each chain has one fermion. 
To make the equations look simpler, let x = e tpi and y = e ip2 , and j and k represent the locations 
of the fermions on chains 1 and 2 respectively; by our conventions j is an odd integer and k even, 
so we have 

<p 12 (j\k) = A\% X iy k + Af x y3x k for j < k 

<f 21 (j\k) = Af 2 x k yi + A 2 2 \y k xi for k < j . 

We then let the Hamiltonian act on the state (pQ) with this (p. Ignoring the boundary conditions 
momentarily, a little algebra shows that this is an eigenstate if 

2p Q (k - l\k) = -ip Q (k + l\k) - iip Q '(k + l\k) - ip Q (k - l\k - 2) + iip Q '(k - l\k - 2) (15) 

and its parity conjugate 

2ip Q '(k + l\k) = -ip Q '(k + l\k + 2) - itp Q (k + l\k + 2) - (p Q ' (k - l\k) + if Q (k - l\k) (16) 

for all even k. Here Q = 12 and Q' = 21, but it is straightforward to show that these constraints 
apply in general when Q' is defined by reversing Qa and Q(a+l) in Q, i.e. Q' = (. . . Q(a— 1), Q{a+ 
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1), Qa, Q(a + 2), . . . ). Then (fT5j) and (fl~6j) must hold whenever Qa < fi and Q(a + 1) > /i, i.e. 
JQa is odd (is on chain 1), and jQ( a +i) is even (is on chain 2). 
The eigenvalue for such an eigenvector is simply ©, i.e. 

E {1>1) =A + x 2 +x - 2 + y 2 + y- 2 = ( x + x" 1 ) 2 + (y + y^ 1 ) 2 . (17) 

Plugging the Bethe ansatz (|13j) for the eigenvector into the constraints (|15|16p gives two equations 
for the four unknowns We find 

p f A n\-( " 2 (* + x ~')(y + i{{x + x^f-{y + y^f)\(A^\ . 

{<) ~ U(* + x- X ? ~ (V + y- 1 ) 2 ) "2(x + x- 1 )^ + y- 1 ) J {AH) { ' 

When (y + y^ 1 ) = ±i(x + x^ 1 ), we have ^(1,1) = 0, and the particles form a Cooper pair. Since 
the matrix has determinant (E^ ^) 2 , it is still possible to solve ()18f) for a Cooper pair, yielding 
A 12 = ±A 12 and j4 21 = j4 2 i = 0- 

This eigenstate must also satisfy the boundary conditions. For periodic untwisted boundary 
conditions, they are simply 

(xy) L = 1, y L A\ 2 = A 2 \, x L Af 2 = Af 2 . 

The first of these is simply the requirement that the total momentum be an integer multiple of 
2tt/L, while the others are the individual momentum quantization conditions in the presence of 
interactions. These equations can be solved to yield all of the eigenstates. The Cooper pairs do 
not occur here, since to satisfy both the periodic boundary conditions and (I18p would require all 
A% = 0. 

The generalization of this computation to arbitrary numbers of particles is now standard, even 
though it is a mere 39 years since the invention of the nested Bethe ansatz. The key ingredient is 
the i?-matrix, which encodes the relations on the necessary to make the state an eigenvector. 
These relations are given in (|14ll5ll6p : one must satisfy the constraints for each pair of fermions 
when they become adjacent. Let P = (. . . Pa, P{a+ 1) . . . ) and P' = (. . . P(a + 1), Pa . . .). Then 
we have 

Ep a ,P{a+l)Ap>,Q = R{pp a ,p P ( a+1 ))Ap : Q (19) 



where 



R(Pa,P b ) 



f-E ab \ 

-8cos(p a ) cos(p b ) 4i(cos 2 (p a ) - cos 2 (p b )) 

4i(cos 2 (p a ) - cos 2 (p b )) -8 cos(p a ) cos(p b ) 

V -E ab J 



(20) 



where E a b is the energy of particles with momenta p a and pb, i.e. 

E a b = 4(cos 2 (p a ) + cos 2 (p b )). 

The i?-matrix is 4 x 4 because there are four possibilities for which chains the fermions at locations 
jQ a and jQ( a +i) occupy. The 2x2 block in the middle applies to the cases where the fermions are 
on different chains, and is the same as that in (fT8|) . The inverse {R{p a -,Vb))~ l = R{Pb-,Va) / {E a b) 2 ■> 
as necessary for consistency. We have followed [10] and written (|19p so that it still is applicable 
when E a b — ► 0, i.e. when p a and pb satisfy the Cooper pair condition ((HJ). 

The R matrix relates coefficients A^, to A^. Repeatedly applying (fT9l) gives all the A® 
in terms of those for a given permutation, which we label as P = 0. For periodic boundary 
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conditions the A® are then related to each other by imposing the boundary conditions as well. 
To derive the resulting Bethe equations in situations like ours where the R matrix is non-diagonal, 
one must do a second Bethe ansatz, the "nested" part of the nested Bethe ansatz [JJ. Here this 
requires one level of nesting, and we omit it here because we will not need it. A vital consistency 
condition in applying (|19|) is that the i?-matrix satisfy the Yang-Baxter equation, which ensures 
that this procedure defines all the A® consistently. This i?-matrix (|20p indeed satisfies the Yang- 
Baxter equation, and for example arises in vertex models [TS] and the t — J model with the SU (2) 
symmetry deformed to SU{2) q [16]. It also arises in the scattering matrix of the kinks in the 
sine-Gordon model at its supersymmetric point 1 = 167r/3 [T7j. We note also that this i?-matrix 
satisfies a special condition called the "free-fermion" condition [18] . which is part of the reason 
behind the miracles in this paper. 

We now return to the boundary conditions of interest in this paper: open ones with a boundary 
magnetic field ([3]). To implement these open boundary conditions, we need to take combinations 
of wavefunctions with =tp a . The Bethe ansatz (|13p is generalized here to 

/ 

P a=l 



Yl A P II (e ipPajQa - 7^ a e- ipPaj ^A . (21) 



The boundary conditions can be satisfied by choosing the coefficients jp a for the momentum 
corresponding to the leftmost and rightmost fermions on each chain. By definition of Q, the 
leftmost fermion is located at j'qi, while the rightmost one is located at jqj. For each Q, the 
form of the boundary conditions depends on which chain the leftmost and rightmost fermions are 
on. We have 



= 
= 



U q Uqi = !,••• !■■■) + V Q (3Qi = -1, ••• | ••• ) if Ql < fi 

\v Q (...\j Q i = o,...) if Ql > /i 

'<p% ..{..., jQf = 2L) + <p% . . | . . . , j Qf = 2L + 2) if Qf > fi 
yj(...,j 0/ = 2L + l|...) ifQ/</i 

Plugging (|2ip into these boundary conditions, we find the same constraint for either condition. 
We have for the left end 

7?i = "I (22) 
for all Q. Likewise, whichever boundary condition for the right end is applicable, we find 

e -i Pp/ 2(2L+l) 7 ^ = _ L (23) 

First let us consider the case where ^ for any choice of a, b, i.e. there are no Cooper 
pairs. Deriving the Bethe equations and the relations among the coefficients now proceeds 
in the same fashion as for periodic boundary conditions. In fact, the computation is virtually 
identical because the i?-matrix is independent of the signs of p a and pb- Thus we can expand out 
the product in the ansatz (|2ip , and apply the same relation (|19p to any of the resulting 2^ terms 
for a given Aq. Consistency then requires that all the jp a for a given Q to be the same. Since 

one of them is fixed to be —1 by the boundary conditions at the left end, this means jp a = — 1 
for all P, a, and Q. We still need to satisfy the boundary conditions (|23p at the right end. Since 
these must hold for all permutations P and Q, the only way this is possible is if 

e ip a 2(2L+l) = 1 / 24 \ 
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for all a. This is precisely the momentum quantization condition for a free fermion 

7T 

Pa = "la 

F 2L + 1 

given in ©! 

Another miracle is that there are multiple states with the same energy: fixing the p a does 
not fix all the Because the boundary conditions are the same for all Q, they do not restrict 
the A^> at all. The relation (|19p can be used to relate the vlp to those for a given permutation. 
Labeling this given permutation by P = 0, (|19p determines the other in terms of A® , but 
it does not fix the latter. The only relation between the different Aq is therefore that required 
by Fermi statistics for the fermions on each chain. Fermi statistics for chain 1 means that when 
both Qa < fi and Qb < /i, we have A® = —A® , where Q 1 differs from Q by swapping Qa and 
Qb, i.e. if Q = (. . . , Qa, . . . , Qb, . . . ) then Q' = (. . . , Qb, . . . , Qa, . . .). Fermi statistics for chain 2 
means that when both Qa > f\ and Qb > $2, A® = —A® as well. For a given set of momenta 
satisfying (|24p . there are therefore (f 1 ^/ 2 ) undetermined coefficients in the Bethe ansatz. All of 
these are perfectly valid eigenstates, and they all have the same energy. We thus recover the 
degeneracy d,E in fjXOj) in the special case where there are no Cooper pairs. 

The natural way of interpreting these results is for the quasiparticles to be the exclusons 
described in section[2j These quasiparticles have a statistical interaction, because each momentum 
p a must be different: the wavefunction vanishes if p a = p for a ^ b, because R(p a ,p a ) = — 1. 
Moreover, even though fermions cannot hop between chains, the wavefunction is not a product of 
separate factors for each chain. A given momenta is not associated with one chain or the other: 
the sum over P includes permutations over all / = f\ + fa momenta. These facts are exactly 
what happens with two species of quasiparticles obeying exclusion statistics. 

We now can see explicitly that we must include Cooper pairs to make the Bethe ansatz 
complete. There are (^) choices of momenta, each with degeneracy (J ) . However, 

;.)©>(/)(/, 

so there are more states than those involving solely exclusons. The derivation that all momenta 
must obey the free-particle quantization condition (f24"j) is valid when E ao ^ for all p a and pb- 
Therefore the missing states must include at least one Cooper pair, a state where E ao = for at 
least one pair (a,b). 

So let us first return to the case where f% = fa = 1, and now let p\ and p2 satisfy the Cooper- 
pair condition ([8|), so that E = 0. As noted above, the requirements (|15ll6p are satisfied here if 
A\% = Af^, and A^f = A%\ = 0. The wavefunction of a single Cooper pair is therefore 



■ A(xi - x-i)(y k - y 4L + 2 " fc ) j < k 
A{x* - x ){y J -y + 3 ) k < j 



where (y + y^ 1 ) = i(x + x _1 ). One can check directly that this yields an eigenstate of Hf + H c 
with eigenvalue for any value of x. Note that the momentum quantization condition (|24p does 
not apply to the momenta in Cooper pairs. 

To show that all the missing states are given by Cooper pairs, we count how many linearly 
independent eigenstates ([25]) yields. It is useful to define the new variable z = (x + = 
i{v + V Choosing 

y-(2£+i) 

A 



(x-x l ){y-y l ) 
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gives ipc as a polynomial in z 2 (it can be expressed in terms of Chebyshev polynomials explicitly 
if desired). It is simple to then see that the highest order of z appearing in these polynomials is 
2 2 ( L_1 ) ) and the lowest is z°. Since z 2 is arbitrary, it follows that the dimension of the space of 
states of a single Cooper pair is L. This agrees with the general degeneracy formula for C = 1 
and no exclusons, N\ = N 2 = 0. This gives all the missing states for fi = f 2 = l: there are 
L(L — 1) two-excluson states, giving the correct total of L 2 . Thus we have the nice result that 
Cooper pairs make up the states missing due to imposing the exclusion statistics. 

Now we consider the general case. It is important to note that we do not single out two 
fermions and call them a Cooper pair. Rather a Cooper pair occurs when two momenta obey the 
condition (|8|). The relation (|19p still is sufficient to satisfy the eigenvector requirements f 1 1 5 II 6 1) : 
(I14j) is satisfied as long as we divide out any E a b from both sides before taking them to zero. 

When Cooper pairs are present, (|19|) means that some of the vanish, and others must 
be equal. For C Cooper pairs, we order the p a so that cos(p2c-i) = icos(p2 C ), for c = 1, . . . C. 
Given a permutation P, define a\ and 02 so that Pa\ = 2c — 1 and Pa2 = 2c. Then for A® not 
to vanish, we find that for all c 

1. (2c — 1) must appear before (2c) in P, i.e. a\ < a,2 

2. if there are no fermions between those with the Cooper pair momenta (i.e. 02 = a\ + 1) 
then these two particles must be on different chains. This means Qa\ < f\ and Qa,2 = 
Qifli + 1) > fi, or vice versa. 

In the latter case, the i?-matrix relation (]19p also requires that where Q' is Q with 

a\ and 02 = o-i + 1 reversed. 

Applying (|19|) means that the free-particle momentum quantization condition (I24p is still 
applicable for all the momenta not part of Cooper pairs, i.e. p a for a > 2C. However the momenta 
in Cooper pairs need not satisfy it. The reason is as follows. Because Pa\ = 2c — 1 must always 
appear before Pa2 = 2c in a P with a non-vanishing , the fermion at position jQ ai must always 
be to the left of a fermion at position jQa 2 ■ Thus no fermion with momentum p2 C can ever be 
the leftmost particle, and no fermion with momentum P2c-i can ever be rightmost. Therefore, 
the boundary condition (f22l) does not apply to any j 2c , and the boundary condition (f23l) does 
not apply to any ^c-i- Since (f24"|) arose from demanding that the two boundary conditions on a 
given momentum be consistent, it does not apply to any Cooper-pair momentum. 

Let us illustrate these conditions with f\ = 2 and /2 = 1, where there are 6 different permu- 
tations for P and three for Q to consider. (We only need consider (123), (132), and (312) for Q 
because the other three are given by applying Fermi statistics to the two fermions on chain 1.) 
With no Cooper pairs, there are three independent A® 23 , as discussed before. For one Cooper 
pair, we have p\ and P2 satisfying ([8]). This means that when P = (213), (231), (321), A^ vanishes 
for all Q. When P = (123), we have a± = 1 and 02 = 2, so this means ^ijJ! = as well- The 
non-zero ones must obey for example ^123 = ^123- Thus for the three A® 23 , one is zero and the 
other two are equal. Applying (|19p then determines all the other A% , and including the boundary 
conditions requires that p% obey the free-fermion condition ([5]). Therefore all the A% are given m 
terms of ^123 : all the degeneracies here arise from the arbitrariness of p\ . Counting the different 
degrees of freedom here is similar to the two-particle case. We replace p\ and P2 with z, so that 
the wavefunction is given in terms of polynomials in z. However, because of the presence of the 
excluson with momentum p^, the lowest-order term in </? 312 is now order z 2 , so there are only 
L — 1 terms in the polynomial. When there are no other exclusons, there are L terms. The reason 
for this excluded volume effect is Fermi statistics. The wavefunction must vanish when j\ = j'2, 
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so neither p 2 nor p\ can be p%, thus reducing the number of possibilities for the momentum. This 
yields a total L{L — 1) states with one Cooper pair (the factor of L coming from the L possible 
values of ps). Likewise, there are 3(g) states with no Cooper pairs, the factor of 3 coming from 
the independent ^4^ 2 3) ■ This is a complete set of states: 

3 f )+L(L-l) 



1 / V 2/ V 3 



To obtain the general degeneracy given in (jlOp . one proceeds in this fashion. The first factor 
comes from the number of independent Aq , while the second comes from the order of the polyno- 
mials in z c describing the momenta in the cth Cooper pair. Let us first discuss the different A® , 
where is our label for (12 .../). Applying Fermi statistics to chain 1 relates all the Q obtained 
by exchanging any of Qa with Qa < /i, while applying it to chain 2 relates any obtained by 
exchanging any of the Qa with Qa > f\. Thus the largest possible number of distinct Aq is (/). 
We showed above that these are all distinct when C = 0. Now consider C = 1. The second condi- 
tion for non-vanishing A® described above means that Aq vanishes for any Q where jq\ and jQ2 
are either both even or both odd. The non-vanishing ones therefore have Ql < f\ and Q2 > f\ 
or vice-versa. There are 2/1/2(7 — 2)! such possibilities. If Q' is Q with Ql and Q2 reversed, 
then Aq = Aq , so the number of distinct possibilities is lowered to /i/2(/ — 2)!. Including the 
effects of Fermi statistics as well means dividing by f±\ and / 2 !. Thus the total number of distinct 
Aq for one Cooper pair is 1) • For general C, applying the second condition C times in the 
analogous fashion gives 

h\ f 2 \ (f-2C)\ (f-2C 



(/i-C)!(/ 2 -C)! h\f 2 \ V/i-C, 

for the number of independent Aq . Thus we indeed recover the first factor in (|10|) . 

The second factor likewise follows from essentially the same argument we gave for small f\ 
and / 2 . The presence of fermions with momenta not obeying the Cooper-pair conditions creates 
an excluded volume effect reducing the number of independent terms in the polynomials in z c . 
The Fermi statistics means that the wavefunction must be symmetric under exchanges of the z c , 
implying that the number of linearly-independent choices of the z c is the binomial 

L - (/ - 2C) 
C 

By writing the grand canonical partition function as at the end of section [21 we saw that including 
the degeneracy (fTUI) successfully counts all the states in the theory. Thus the Bethe ansatz is 
complete. 



4 The extended (super)symmetry algebra 

The extensive degeneracies in the spectrum given by (fTUj) imply that our model must possess 
an extended symmetry algebra. This is obvious in the quasiparticle basis of states in terms 
of exclusons and Cooper pairs. A symmetry operator creates or annihilates a Cooper pair, or 
changes the type of a excluson, preserving the momentum. In this section we describe this 
symmetry algebra, and how it acts on the orginial fermionic states. The symmetry turns out to 
be an extension of supersymmetry, and is infinite-dimensional as L — > 00. As we will show, it is 
quite reminiscent of the Yangian symmetry of the Haldane-Shastry model. 
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The basic supersymmetry operators act on the exclusons: Q changes type 1 to type 2, and 
the inverse. Changing the type of an excluson requires flipping a fermion from one chain to the 
other. Thus Q hops a particle on chain 1 to chain 2, with doing the reverse. These operators 
are non-local: they depend on the locations of the other particles. Let 

k 

3=1 

We then define 

L-l 

q = 4<* + E ( e "* 2fc - W2 4 + ^ W2 4 +2 ) caw. (26) 

k=l 

Both Q and annihilate the Cooper pairs. 

This supersymmetry charge commutes with H = Hf + H c if Hf, is the boundary magnetic 
field in ([3]). This can be verified directly, but it is more illuminating to work out the entire 
supersymmetry algebra. Even though Q is written in terms of fermion bilinears, it is effectively 
fermionic because of the non-local string of operators from the exponential. One finds for its 
anticommutators 

Q 2 = 0, (Q^ 2 = 0, H = {Q,Q^}. (27) 

It then follows from (J27|) that [Q,H] = [Q^,H] = 0. 

One can derive a number of interesting properties of the states directly from the (unextended) 
supersymmetry algebra (127p . Since we have already derived them (and many more) by using the 
Bethe ansatz, we will just state these properties here; see [19] for detailed explanations of the 
methods. All eigenvalues E of the Hamiltonian obey E > 0. States with E > form doublets 
under the supersymmetry algebra, while the E = ground states are annihilated by both Q and 
Q>. The Witten index provides a lower bound on the number of E = ground states [20j . Since 
the supersymmetry generators are fermionic, for purposes of computing the Witten index, each 
fermion on chain 1 has charge —1/2, while each on chain 2 has charge +1/2. Then 

W = e i(/2_/l)7r/2 = 2 L , 

states 

so there are at least 2 L ground states for all possible values of /i and fi for a given size L. We 
showed above that a state with E = must be comprised solely of Cooper pairs. These occur 
whenever /i = /2, and from (110p we see that the number of such states for a fixed /i and f2 is 
because C = fi = /a here. Thus there are exactly 2 L ground states, and the lower bound from 
the Witten index is the exact number here. One can find a Hamiltonian with twisted periodic 
boundary conditions obeying (|27p . The Witten index is the same, so the Cooper pairs remain 
for the twisted model, but the other degeneracies no longer occur. 

To understand the full degeneracies, supersymmetry alone is not enough. Q does not change 
the total fermion number, whereas annihilating or creating a Cooper pair will annihilate or create 
fermions on both chains. Starting from a state with all particles on the upper chain, acting with 
Q flips one particle to the lower chain. Since Q is nilpotent the action cannot be repeated. 
Moreover, To obtain the full degeneracies (/) of the exclusons one would need / generalized 
supercharges that can act independently. We have found that a hierarchy of operators Qp, Q~ , 
Hp, with p = 0, 1, . . ., can be constructed with the following algebraic properties 

{Q+,Q+} = 0, {Q-,Q~} = 0, {Q;,Q-} = H p+q , (28) 
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implying that [Qp, H q ] = for all p, q. We refer to this algebra as a 'super-Yangian', by analogy 
with the Yangian appearing in the Haldane-Shastry chain [12j. 

The leading terms in the hierarchy are Qq = Q, Qq = , Hq = Hf + Hj, + H c . For the 
operator Qf we choose a form where a particle at site % in the upper chain hops to sites i + j, 
j = 1, 3, ... in the lower chain with amplitudes containing an alternating sign (— 1)0'— W 2 ^ 

Qt = E e ~ ia2k7T/2 ( E (-l)°'- 1)/2 ^ j7r/2 4 + i + , J c 2k+ i , (29) 

k=0 \j=l,3,...,27V-2fc-l / 

with ct2k as given above and 

02k,j = 2(n 2fc+2 + "-2fc+4 + • • • + n 2 k+j-i) ~ (n 2 k+3 + n 2k+5 + • • • + n 2k+j ) ■ 

On a 1-particle Bethe state with x = e tp these amplitudes combine as 

aT 1 - x" 3 + x~ 5 (a; + x -1 ) -1 . 

As a consequence, we find that on a general 1-particle Bethe state the operators H p act diagonally 
with eigenvalue 

E p [x] = (x + x- 1 ) 2 ^ . (30) 

On a general /-particle state we find that H\ = for / even while H\ = 1 on all states with / 
odd. 

We have investigated the explicit form of with p > 2 and H q with q > 3 for small system 
sizes, confirming that operators satisfying the full algebra eq. (j28j) indeed exist. Restricting to 
2-particle states, there is some freedom to choose the ; what seems to be a canonical choice 
leads to eigenvalues 

E p [x, y] = (x + x- 1 ) 2 - 2 ? + (-lyfo + y- 1 ) 2 - 2 ? (31) 

on states with two exclusons. In addition, the H p have an L-fold degenerate eigenvalue E p = 
when acting on states with a single Cooper pair. States with more particles show a similar pattern; 
for / particles with C Cooper pairs the eigenvalues E p reduce to those of / — 2C exclusons. The 
operators can be chosen such that on a state with exclusons at {x a }, a = 1, . . . , /, the eigenvalues 
of Hp with p even take the form 

/ 

Ep[{x a }] =^2(x a + x^) 2 - 2p for p even. (32) 

a=l 

The algebraic structure in eq. (j28H is highly reminiscent of the Yangian Y(sl 2 ) which features 
as the symmetry algebra of various integrable models [21]. In the Haldane-Shastry spin-chain, 
the S'L r (2)-spin symmetry, with generators Qq with A = ±,3, extends to a Yangian generated 
by Qp, P = 0, 1, . . ., with all Qp commuting with conserved quantities H q [12]. This algebraic 
structure beautifully reflects the underlying physical picture of 'spinons', with the Yangian gen- 
erators performing spin-flips on a multi-spinon state with all spins polarized. In our model 
here, the situation is similar, with the Q^ sweeping out a full multiplet of eigenstates from a 
reference state with all particles 'polarized' on the same chain. Comparing with the Haldane- 
Shastry model, the role of the S'f7(2)-spin symmetry is taken over by supersymmetry, bringing 
the charges Q^ and the conserved quantities H q together in one algebraic structure. The actual 
form of the first nontrivial 'sYangian' generators Qf is reminiscent of the form of the (bosonic) 
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Yangian generators (sometimes referred to as the 'logarithmic Yangian') in integrable quantum 
field theories of massive particles and WZW models of conformal field theory [22] . 

The results of this section apply to open boundary conditions with the boundary magnetic 
field ([3]). An interesting open question is if this symmetry algebra can be deformed to apply to 
other boundary conditions, where the model is still solvable but there no longer exist the large 
degeneracies (jlOp . If such a symmetry algebra does persist, it is then spectrum-generating like 
the Yangian. 

We remark that in addition to all symmetries already mentioned, the model (with open 
boundary conditions) admits a discrete symmetry reminiscent of that of the Hubbard model 
at half filling. It is implemented by interchanging particles and holes on the upper chain only, 
leading to 

(/i,/ 2 )^(£-/l,/2), E ~2L-1-E , (33) 

and in particular 

N 2 <-> C . (34) 

The latter relation shows that the particle-hole transformation maps the operators (which 
each flip one particle from upper to the lower chain) into operators creating Cooper pairs. 

5 Conclusion 

We have discussed in depth a one-dimensional model with local interactions whose excitations 
are Cooper pairs and quasiparticles obeying exclusion statistics. We used the Bethe ansatz to 
derive these results, and found that various technical complications ended up being crucial to 
describing some very interesting physics. 

Despite the presence of Cooper pairs, our system is not really a superconductor, because it 
is gapless. However, we believe it is possible to modify the Hamiltonian to give the exclusons 
a gap while preserving the Cooper pairs. The reason is the supersymmetry. If we deform the 
couplings without breaking the supersymmetry, the Witten index will not change [20], so the 
2 L ground states with E = will remain. Any supersymmetry-preserving deformation which 
gaps the exclusons will therefore leave the Cooper pairs. The next step would be to break the 
supersymmetry slightly so that Cooper pairs can have a kinetic energy below the gap. 

Our model has some intriguing similarities to a solvable one-dimensional model of a super- 
conductor. The "Russian-doll" model is an effective model of Cooper pairs, where time-reveral 
symmetry is broken by including a complex phase in the Hamiltonian [23J. This is very reminis- 
cent of how time-reversal symmetry is broken in our model, and it would be interesting to find a 
deformation of our model whose Cooper pairs are described by the Russian-doll model. 

We believe our model will prove useful in understanding exclusion statistics as well as su- 
perconductivity, since it can be treated with well-studied Bethe-ansatz methods. It should be 
possible to derive many more interesting properties for it. We are also hopeful that it may shed 
light on anyonic superconductivity. We find it amazing that 75 years after the invention of the 
Bethe ansatz, it is still being used to derive and explore fascinating properties of one-dimensional 
quantum systems. 

We are very grateful to Fabian Essler for his patience. This work was supported in part by 
the foundations FOM and NWO of the Netherlands, and by NSF grant DMR-0412956. 
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